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The light output of the 324 CsI(Tl) scintillators of INDRA has been measured over 
a large range both in energy: 1 — 80 AMeV and in atomic number of incident ions: 
Z = \ — 60. An analytical expression for the nonlinear total light response as a 
function of the energy and the identity of the ion is proposed. It depends on four 
parameters. For three of them, connected to CsI(Tl) intrinsic characteristics, recom- 
mended fixed values are proposed. They are issued from the comparative study of 
the forward scintillators of INDRA, based on intermediate mass fragment data. The 
fourth one, related to light collection and to the gain of the associated photomulti- 
plier, is particular and may be accurately obtained afterwards, from light charged 
particle data. Two applications are presented: fragment identification in telescopes 
using a CsI(Tl) crystal as residual energy detector and the scintillator energy cal- 
ibration. The results are successfully confronted to heavy fragment experimental 
data. 

Key words: PACS number: 29.40.Mc, 32.50.-Fd 

(light response of CsI(Tl) to heavy ions, quenching, delta rays) 



1 Introduction 



INDRA is a Att axially symmetrical array for the detection of light and heavy 
charged nuclear reaction products [1,2] covering huge dynamic ranges, both in 
energy (from ^1 MeV to the maximum available energy ^6 GeV at GANIL) 
and in atomic number (from proton to uranium). It has a high granularity 
and a shell structure, consisting of several detection layers. For the last layer, 
which should stop all particles and fragments produced using GANIL beams, 
thallium-activated caesium iodide scintillators (CsI(Tl)) coupled to photomul- 
tiplier tubes were chosen. 

The standard calibration procedure for CsI(Tl) is to find a function L depend- 
ing on the energy, but also on the identity of the particle, which describes 
reasonably well the induced scintillation Qq. The parameters of this function 
are determined by a global fit procedure which simultaneousely compares the 
calculated scintillator response to the experimental one for all particles and 
fragments of well known energies. Afterwards, the unknown energies of the re- 
action products detected in physical runs arc found by means of the calibration 
function and the related parameters, starting from the measured associated 
light outputs. Once this is accomplished, one may get a reference map for the 

1 Corresponding author. Tel 33 1 69157148; fax 33 1 69154507; e-mail borderie® 
ipno.in2p3.fr 

2 present address: DRFC/STEP, CEA/Cadarache, F-13018 Saint-Paul-lez-Durance 
Cedex, Prance. 
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identification of the reaction products in a two-dimensional plot showing the 
energy deposited in the preceding detection layer, for every fragment punching 
through, versus the total light from the CsI(Tl) crystal. This is the second aim 
of our work. 

For the forward angles of INDRA (3° < ^ < 45°, rings 2-9), the detection 
layers which precede the scintillators consists of gas ionisation chambers (ICs) 
and 300 /xm silicon detectors. The Si detectors have allowed an accurate de- 
termination of the residual energy - as presented in subsection 2.2 - for the 
whole range of fragments passing through and stopped in the scintillators. 
Thus, there is a tremendous set of data which has facilitated a detailed study 
of the Csl(Tl) crystal light response, having led to a convenient expression 
which describes the total scintillation. 

Both components of the procedure, measured scintillation and its modelling 
function, have implied special processing. In the case of INDRA, neither the 
amplitude of the CsI(Tl) associated signals nor the integral of these signals are 
measured. Instead, fractions of the total signal are integrated into two time 
gates, allowing particle identification. Both integrals have rather complex de- 
pendences versus the incident particle energy as compared to the whole inte- 
gral. It is possible to accurately find the latter quantity by software, starting 
from the two measured signal fractions, as shown in sections 2 and 3. 

The exact expression of the total light output of a CsI(Tl) crystal, as predicted 
by the recombination and nuclear quenching model (RNQM) [3], implies a nu- 
merical integration over the energy and this fact is prohibitive for application 
purposes. Under suitable approximations, the integration may be however an- 
alytically performed and a very easily to handle light response expression is 
deduced. Starting from intermediate mass fragments of known energies, the 
3 or 4 involved parameters are determined. Except for the gain parameter, 
proper to each scintillator crystal and corresponding electronic chain, the val- 
ues of the other associated fit parameters are fixed, as being connected to 
intrinsic CsI(Tl) crystal properties. Procedures for fragment identification in 
a AE — Qq telescope-type map, with the scintillator as residual energy de- 
tector, as well as for the energy calibration of the latter one - in the whole 
dynamic range - are developed and critically analyzed in section 4, containing 
the RNQM apphcations. 

For the backward angles of INDRA (45° < 6 < 176°, rings 10 - 17), the 
scintillators are preceded only by ionization chambers. The calibration of the 
scintillators leans on the above mentioned light response expression and the 
CsI(Tl) characteristic parameter values, found at forward angles; the indi- 
vidual gain parameter is determined by means of light charged particle and 
eventually light fragment data. The calibration so found allows to safely ex- 
trapolate the charge identification in a AE — Qq map to regions where no ridge 
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lines arc visible because of very low statistics, improving then both charge and 
energy determination for heavy fragments {Z > 15) detected beyond 45° with 
INDRA. Details are given in section 4. 

Our findings are summarized in section 5. 

Notation and values of physical constants and variables used in this paper. 
See also those in the preceding paper [3]. 
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Symbol* Definition 



Units or Value 



Experimental light output and related variables 



t 

i{t) 

Qf 

V 
Qs 



Qfs 
Qo 

To 
T 

F 



ai 



02 



03 



04 



fgeom 
fpMT 



time 



a.u. s 



signal at the last dynode of the CsI(Tl) 
photomultipher (PMT) 

integral charge of the fast component of a.u. 
the signal 

decay time constant of the fast component s 

integral charge of the slow component of a.u. 
the signal 

decay time constant of the slow s 
component 

integral charge corresponding to the whole a.u. 

signal 



measured signal at the output of the PMT 
anodic circuit 



a.u. s 



approximate total integrated charge oc ex- a.u. 
perimental light output 

decay time constant 
lower value of the decay time constant 
rise time constant at output of the PMT 



s 
s 
s 

a.u. 



experimental charge integrated in the 
"fast" gate 

experimental charge integrated in the a.u. 
"slow" gate 

Calculated light output and related variables 

constants in the approximative expres- a.u. 
sions of Se,n 

gain fit parameter in the friendly analyti- a.u. 
cal expression of L 

quenching fit parameter in the friendly an- a.u. 

alytical expression of L 

fit parameter in the friendly analytical MeV 
expression of L 

fractional energy loss transferred to a S - 
ray, a fit parameter in the friendly analyt- 
ical expression of L 

light collection factor a.u. 

PMT gain factor a.u. 
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*Most of the notations of the original references have been kept. 



2 The thallium-activated caesium iodide scintillators 

2.1 CsI(Tl) crystals of INDRA and associated electronics 

There are 324 CsI(Tl) crystals in all, with thicknesses ranging between 138 
mm and 50 mm from forward to backward angles [1] . All the crystals belong- 
ing to one ring have identical shape and size. Because the light output of the 
CsI(Tl) crystals critically depends on the temperature [4,5], a water cooling 
system stabilizes at 20° C the temperature in the mechanical supports of the 
CsI(Tl). The CsI(Tl) crystals are coupled to photomultipliers tubes (PMT) 
[1]. The use of PMTs provides lower energy thresholds for mass identification 
as compared to those obtained with photodiodes [6]. The stability control of 
the scintillators is ensured by optically connecting them to a laser system [1], 
which makes use of a nitrogen laser [7] , emitting in the ultraviolet (UV) band, 
and CsI(Tl) hght distributors. The CsI(Tl) PMT signals are fed in 24 input 
VXI bus modules containing the processing functions. Each channel comprises 
a constant fraction discriminator, two integrators for "fast" and "slow" parts 
with accompanying delay and gate generators. The analog to digital conversion 
is performed by two multiplexed 12 bit converters. Exhaustive descriptions of 
the CsI(Tl), PMT and associated electronics, as well as of the data acquisition 
and triggering system are given in [1,2]. We are doing here only a short pre- 
sentation, stressing those details which are necessary to make comprehensible 
the procedures concerning the CsI(Tl) energy calibration and the fragment 
identification in the IC - CsI(Tl) maps of the backward rings. For energy cal- 
ibration purpose, rings 10 - 17 were each equipped with a single two-element 
telescope (80 jivn and 2 mm thick sihcon detectors) - which will be referred to 
as the calibration telescope (CT). The CT covers part of one of the CsI(Tl) 
crystals in each backward ring. 

2.2 Calculation of the deposited energy into the scintillators 

For rings 2-9, the calculation of the energy deposited into the CsI(Tl) crystal 
to which a given hght output corresponds - E'o - is based on the energy lost in 
the preceding layer of a detection cell of INDRA, using the nominal thickness of 
each silicon detector (si SOOyum). The silicon detectors were carefully calibrated 
taking into account the pulse height defect, with an absolute accuracy of 
2 — 3%. [8]. However, the relative accuracy - between different ions up to Xe 
and for different energies up to 80 AMeV - is within 1%. For ions as light 
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as Boron, the consequent relative accuracy for the residual energy in CsI(Tl) 
crystal is 1.3% at 50 AMeV, 1.4% at 25 AMeV and 3% at 5 AMeV. For ions 
as heavy as Xenon, the consequent relative accuracy for the residual energy 
in CsI(Tl) crystal is 1.8% at 50 AMeV, 2.5% at 25 AMeV and 10% at 5 
AMeV. One has to note that the relative accuracy of the total energy of the 
ion, deposited in both Si and CsI(Tl) detectors never exceeds 2 — 3%. Starting 
from the energy deposited in the 300/um Si detector, and by using the stopping 
power tables of Hubert et al. [9] above 2.5 AMeV and the renormalized variant 
of those of Northcliffe and Schilhng [10] under 2.5 AMeV, the residual energy 
deposited in the scintillator placed behind it is calculated. For light charged 
particles and light fragments (Z < 4), the mass is identified, while for heavier 
fragments, an hypothesis is necessary for the mass. 

For rings 10 - 17, the CTs play the essential role in estimation of the energy 
deposited into the scintillators, as providing reference energy spectra. Instead 
of directly looking for the correspondence: total light output - deposited en- 
ergy into the scintillator, event by event, the incremented spectrum (for each 
reaction product) is compared to the associated reference one, to which it is 
stretched. The calibration procedure for the backward rings is presented in 
subsection 4.4, followed by a description of a rapid fragment identification 
recipe making use of the same mentioned expression. 

All along this paper, the Xe + Sn system at 32 and 50 AMeV incident ener- 
gies is used and emitting reaction products over a large range, in energy and 
atomic number, have been selected. The experimental data were taped only for 
multiplicities higher than or equal to 4 (multiplicity > 4). Thus, the reaction 
products in the present application originate only in exit channels involving 
a non negligible transfer of kinetic energy into internal degrees of freedom. 
As a consequence, the emitted primary fragments have excitation energies 
(greater than around 1.5 AMeV) leading mainly to neutron evaporation. The 
secondary fragments (after evaporation) populate the "attractor" line in the 
map of nuclides [11] rather than the "stability" line. The corresponding mass 
formula [11] will be consequently employed when the isotopic mass was not 
determined. In most of the cases, one CsI(Tl) detector will be used (module 2 
of ring 3) to illustrate the described procedure. The nominal thickness of the 
preceding silicon detector is 304//m. 



2.3 The shape of the signal 

The light emitted by a CsI(Tl) crystal hit by a charged reaction product 
has a rise time negligible [12-14] as compared to the decay time which is in 
the microsecond range. The rise time is related to the transfer of the energy 
deposited by the particle to the optical level involved into the scintillation. 
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while the decay time constant concerns the hght emission. 

Traditionally, the decaying part of a CsI(Tl) scintillation is described by one 
[15,16] or two [13] exponentials of short decay-constant (fa 1/is) or, more often, 
by one short decay-constant exponential (~ l/xs) and one long decay-constant 
exponential (7 /is [17]). The short decay-constant depends on the identity of 
the particle, while the long one is considered to be the same for all particles. 

Proton induced signals recorded by means of flash ADCs up to 20 fis [18] 
have shown decaying parts which are curved in a semilogarithmic scale, with 
a steep descent in the first 3-4 /is, which includes the larger part of the integral 
of the signal. The lower is the incident energy, the steeper is this part of the 
signal and the higher is its weight in the signal integral. For very low energy 
(a few MeV), the shape of the decaying signal is a straight line in a semiloga- 
rithmic representation versus time, i.e. it shows an exponential decaying time 
dependence. At higher energy, several exponentials would be necessary for a 
perfect description of the shape of the whole decaying curve. As in Ref. [17], 
the authors of Ref. [18] have kept only two. The first one - the "fast compo- 
nent" - for the dominant steep descent part, has a short decay-constant (0.5 
IIS - 1 lis), with a strong dependence on the atomic number Z, mass number 
A and incident energy Eq of the particle. The second one - the "slow compo- 
nent" - has a long decay-constant (5 us [18]) nearly independent of the type 
of particle. 

Pulse shape analysis allows particle discrimination up to Z=5 for the Csl(Tl) 
scintillators of INDRA. At the highest energies in Ref. [18], the fast component 
covers at least fs 65 % of the integral of the signal for hydrogen isotopes at 20 
- 40 AMeV, 75 % for helium ones at 30 AMeV, ^ 85% for hght fragments 
(Z=3-6) at 15 - 25 AMeV [18] and at least 95 % for heavy fragments, as shown 
in Fig. 1 for Si at 8 AMeV and Kr at 50 AMeV (present work). The above 
mentioned weights become even more important when the incident energy Eq 
decreases, i.e. when the average specific electronic stopping power Eq/R{Eq) 
(for which a good estimate, within a mutiplicative factor, is AZ'^/Eq) is high 
enough. R{Eq) is the corresponding particle range, and the estimate of the 
average specific electronic stopping power is derived from the approximation 
of Bethe-Bloch formula —{dE/dx)e oc AZ'^/E. From the above considerations, 
one may assume in the latter case that only one decay-constant To(£'o,A,Z) is 
involved for each event. Its value will be close to that of the fast component. 

A signal at the crystal PMT output, described as in Ref. [18] by two expo- 
nential functions associated to the fast and slow component respectively: 

i{t) = — e T H e ^B, (1) 
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Fig. 1. The shape of the measured hght signals induced by 8 AMeV Si and 50 
AMeV ^^ifr ions in two different CsI(Tl) crystals of INDRA - solid symbols - is 
well described by one exponential decaying curve provided by Eq. (3). 

has the total light response corresponding charge — Qi + Qs got by inte- 
grating i{t) over time between and oo. 

In one exponential approximation, the same signal would be: 

i(t) = ^e-i, (2) 



where Qo approximates the integral of the signal, Qfs- 



3 Reconstruction of the total light output 



In the case of INDRA, only parts F, S of the signal are integrated in the time 
gates - 400 ns and 1600 - 3100 ns, respectively. Let us make the following 
exercise: consider the expressions of F and S provided by the two exponential 
formula of the signal (1) on one hand, and by the one exponential formula 
(2) of the signal, on the other hand. The values of the expressions of F and 
S found in both cases have to be equal. From these equalities, one can derive 
To and the ratio Qo/Qts for the data in Ref. [18], but making use of the time 
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Fig. 2. Qq is the approximate integral of the light signal when one decaying ex- 
ponential shape, of decay-constant tq (Eq. 2), is assumed. Qfs is the true one, 
corresponding to two decaying exponential shape (Eq. (1)). Their ratio increases 
with: a) the estimate of the average specific electronic stopping power c>c AZ'^/Eq 
(Eq in MeV), b) the inverse ratio of tq; romin is a normalization constant. Processed 
data from Ref. [18] . 

gates of INDRA. The results of this estimation are plotted in Fig 2a) against 
the estimate AZ'^/Eq of the average specific electronic stopping power. This 
plot has a predictive character: the maximum error done in the integral of the 
signal estimation would be of about 10% in the case of the most energetic 
protons but much lower for the charged reaction products with Z > 1. Fig. 
2b) shows the ratio Qo/Qfs versus the reciprocal of the decay-constant value 
To of the ion in question; the normalization constant romm is in fact the lower 
measured value for ions at the lowest energies (~ 1 AMeV) and hence the 
highest stopping powers. This kind of plot could eventually be used in order 
to correct Qq. In any case, as long as AZ'^/Eq > 0.4 (e.g. Ar ions with Eq/A 
< 810 MeV/nucleon), Qq estimates Qis within 2%. This is the case for most 
of our data. 



In view of the above argument, we shall suppose in the following that the 
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current at the last dynode of the PMT associated to a CsI(Tl) scintillator 
varies exponentially in time as in Eq. (2), where Qq is the total charge. This 
current, injected in the anodic circuit of the PMT, leads to the measured 
current at the output of the PMT base Ves(^); which may be expressed [19] 
by means of the equation: 

WsW = ^(e-^-e-^). (3) 

To - T V / 



The total light output is proportional to the total charge Qo; r and Tq are 
the rise time and decay time constants respectively, r has been measured for 
the bases of all PMT (60 ns for rings 11-16 and 20 ns for other rings). The 
shape of the signal given by equation (3) is shown in Fig. 1. By integrating 
expression (3) within the gates mentioned above, one obtains: 



Qo 



To - r 
QqTq 

To -T 



e ) - r (l - e-*i/^) 



ro(l-e-*i/-°) - 

g— t2/To Q—ts/ro 



(4) 
(5) 



with ti=390 ns, t2=1590 ns and ^3=3090 ns (the integration gates have un- 
dergone a diminution of 10 ns according to the delay of the signal in the VXI 
cards). The reason why it was necessary to approximate the decaying shape 
of the signal by a single exponential function of time - eq. (3) - is that in 
INDRA experiments, only F and S (the channels of the charge convertors 
corresponding to the two gates) are measured and not the integral of the sig- 
nal. Under this assumption and by means of the two eqs. (4) and (5), it is 
possible to find the two unknown quantities tq and Qo and therefore, within 
a multiplicative constant, the total experimental light output. The measured 
resolutions (FWHM) of the CsI(Tl) crystals for Si of 7.86 AMeV and a parti- 
cles of 21 AMeV are of ^ 3% for F and ~ 4% for S. This leads to an accuracy 
of f» 1.3% for Qo above 10 AMeV. Below a total energy of 10 MeV the ac- 
curacy progressively goes down to 3 - 4%, which roughly corresponds to the 
measured resolution on the total light for a particles of 5 MeV. For protons of 
21 MeV, the resolutions (FWHM) are: ^ 4% for F and ^ 7% for S, leading to 
an accuracy of ^ 2% for Qo above 10 MeV. Below this energy, the accuracy 
progressively goes down to 6 - 7%. 



4 Approximate formula from RNQM model 
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4-1 Analytical integration 



In practical situations, an analytical integration of the total light output issued 
from the RNQM [3] would be more suited. This is possible starting from the 
first order approximation for total light output formula (expression (18) in the 
previous paper [3]): 



L = aG 

-Bp 

+ 



1 



dE 



1 + a^S^{E) + anS,{E) 1 + S^{E)/S,{E) 



1 - T{E) 



X 



dE 



Es 
Eo 



-It 

Es 



1 + a^S^{E) + aj,S,{E) 1 + S^{E)/S,{E) 

T{E) dE 
+ S^{E)/S,{E) 



(6) 



if suitable approximations are made for the stopping powers, the concentration 
Nji{E) of the defects created by the incident fragment and the fractional energy 
loss J^{E) deposited outside the primary column by the generated 6 - rays. 
All these quantities are discussed in the preceding paper [3]. 

i) For the specific electronic stopping power formula of Bethe-Bloch. the usual 
approximation: (dE/ dx)c(E) = CcAZ^/E, reasonable above a few AMeV, 
may be used; here is a constant including the logarithmic term in the 
Bethe-Bloch formula, much more slowly varying than 1/E. 

ii) The created defect concentration Nn{E) estimated by A^Ruth(-E'), is well ap- 
proximated by neglecting the second term (oc E"'^) in equation (4) of the 
preceding paper [3] A^n oc AZ^/E. 

iii) The specific nuclear stopping power, (dE/dx)n(E), may also be roughly 
estimated by an AZ^/E behaviour, as shown in Fig. 2b) of the preceding 
paper [3]: { dE / dx)^{E) — C^AZ'^/E with constant. In this way, the 
factor (1 + Sti{E) / S-e{E)) in the denominator of all the terms in expression 
(6) becomes a constant: 1 + Cn/Cc, to be included in the multiplicative 
parameter aQ that will be called ai, and the nuclear and recombination 
quenching terms to the denominator of the first two terms (concerning the 
primary column) in the same expression may be summed and replaced by 
only one: a2AZ'^/E. 

iv) By keeping the zero and first order terms in the Taylor expansion around /9| 
of the logarithmic term of the fractional energy carried by the 6 - rays (see 
Eq. (7) of the preceding paper [3]), one may get an approximate expression 
of ^(/52): 
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Table 1 

Fit parameters oi, 02, as, 04. The errors on the parameters (one unit on the last 
digit) are only statistical. The analytically integrated expression (9) of the total 
light output has been used: a) values obtained by means of data available in the 
whole ion range, for a forward module {9 = 4.5°) of INDRA; b) values averaged 
over 8 modules placed on the forward rings {9 < 45°) obtained only by means of 
intermediate mass fragments and light charged particles. These recommended values 
of a2, as, 04 to be used in Eq. (9) are suitable for all INDRA CsI(Tl) crystals. 







a) 


b) 




[a.u.] 


19.5 


variable 


02 


[a.u.] 


0.71 


0.25 


03 


[MeY/u] 


3.8 


3.1 (1.0)* 


04 




0.26 


0.27 



* see the text for explanation 



- 1 



2 In 



-/3|) + f-l' 



(7) 



With the above items i)-iv) assumptions, the first order approximation for- 
mula (6) of the total light output depends on three parameters only and may 
be analytically integrated. The quality of the fragment loci reproduction in a 
AEsi — Qq map will be shown in the next subsection. 



4-2 A friendly analytical formula for the total light output 



The alternative to the item iv) approximation of the fractional energy carried 
by the 5 - rays is to consider it as a step function of energy: 
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200 400 600 800 1000 1200 1400 

Eo (MeV) 



Fig. 3. Total light output against the energy for different ions (AZ = 5): the symbols 
are experimental data from the system Xe + Sn at 32 and 50 AMeV; the curves are 
calculations done with the simple light output formula Eq. (9). 

0, E/A < as , . 

04, E/A > as ' ^ ' 



HE) 



where 03 is the energy per nucleon threshold for the 5 - ray production and 04 
will be a fit parameter too. The advantage is that the first order approximation 
of total light output expression (6) becomes a simple, easily to handle one: 



ai { Eq 



AZ\ 
a2^z — In 



E^ 







a2 



(9) 



{Es = A y. as), very suitable for energy calibration purposes. The fit param- 
eter values are given in column a) of Table 1, and the quality of the fit is 
shown in Figs. 3,4 (solid lines) and 5. Even if the total light outputs are no 
more as nicely reproduced as by exact calculations [3], especially for high spe- 
cific electronic stopping power values, the description of the reaction product 
identification in the two-dimensional plot (Fig. 4) and the deviations of the 
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Light output (a.u.) 



Fig. 4. A AEsi — Qo map (module 2, ring 3 of INDRA) from the sytem Xe + Sn at 
32 and 50 AMeV. The symbols are experimental data. For the calculated total light 
output, the approximations i) - iv) from subsection 4.1 have allowed an analytical 
integration, by making use of two approximate expressions of J-": Eq. (7) - dashed 
lines - and the step function Eq. (8) - solid lines; the dotted lines are obtained if 
the 5 - rays are completely neglected {T = 0) (subsection 4.5). 

calculated energies relative to the true energies (Fig. 5) remain comparable to 
the exact calculation case. More precisely, the heavy fragment identification 
in a AE — Qq map is possible with a resolution of one unit charge around 
Z=40. The corresponding energy deviations may locally reach up to 15% - 
20%, but globally there are inside ~ 6%. About 3% of accuracy are lost as 
compared to the exact calculations. Note in Fig. 4 that the step function ap- 
proximation for (solid lines) does not worsen the result as compared to the 
physical approximation (7) (dashed lines). In fact, it is approximation i) for 
the Bethe-Bloch formula which is distorting - via the fit procedure - the shape 
of the light response for the heaviest fragments (Fig. 3). 

For the energetic light charged particles, there are discrepancies at high energy, 
whose origin is the slight underestimation of total light output (see sect 2.3). 
For energy calibration purpose, the recipe to ameliorate the situation was to 
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Fig. 5. Deviations (in %) of the energy values, determined with the simple light 
output formula Eq. (9), from the true energy values for several products of the 
reactions Xe + Sn at 32 and 50 AMeV: a) deviations plotted against the product 
energy per nucleon: symbols; the regions between curves show the accuracy of the 
true energy per nucleon (see subsection 2.2 for details); b) deviation histograms. 



use a different gain parameter oi for low average specific electronic stopping 

power {{AZ'^/Eq < 0.4), or directly for protons, which constitute most of 
the data for which the experimentally determined Qq differs from the real 
experimental light output by ~ 2 — 10% (Fig. 2a)). The alternative could be a 
previous correction of Qq concerning these data, based on the Tomin/7o values, 
as suggested by Fig. 2b). 

An interesting point is that the gain parameter values got by using only light 
fragments, or even light charged particles are very near the values obtained 
by using a large ion data range up to Z—A5, if the other parameters, provided 
by the whole range of ions, are kept constant at the values given in Table 
la). Such a conclusion clearly appears in Fig. 6: 04 values are found nearly 
independent of the upper limit in Z (Zmax) considered, when the last three 
parameters are fixed . If only a^, are kept constant at the mentioned values, 
oi, a2 remain practically the same as long as Z^ax > 14. Otherwise, as the light 
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Fig. 6. Variation of the parameter ai as a function of the upper Umit of Z considered 
for fit; 02, a^, 04 were kept fixed with values from Table la). The error bars are 
smaller than the symbols because only statistical errors on the parameters have 
been considered. 



response induced by intermediate mass fragments or light charged particles is 
less non-linear than that corresponding to heavy ions, the parameter 02 has 
the tendency to diminish and, consequently, ai too, for the same quality of 
the description. 
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A search of the fit parameter values was performed for one module of each of 
the 8 forward rings equipped with preceding silicon detectors, by using as input 
data the light outputs induced by intermediate mass fragments, accessible in 
the covered angular domain {9 < 45°) for the studied system. It has led to very 
similar values of the parameters 03, 04, close to those provided by large Z and 
Eq scale data. The quenching parameter 02 is related to the average stopping 
power of the reaction product, i.e. to its identity and energy, but also to the 
activator concentration of the crystal, as we shall see in the next subsection. It 
may vary from one to another ring, but not dramatically. The corresponding 
averaged values of these 3 parameters - presented in column b) of Table 1 - are 
the recommended values for the applications when formula (9) is used. 03 was 
decreased at the lower detection threshold (the value in parenthesis) in order to 
avoid the discontinuity in the energy spectra induced by the "non-derivability" 
at E — Es of the light output expression as a function of the energy. By keeping 
02, 03, 04 fixed, the remaining gain parameter may be accurately determined 
as a free parameter by using simply light charged particles. The results are 
very similar to those shown in Figs. 3,4 with solid lines and in Fig. 5. 
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4-3 Comparative study of the CsI(Tl) crystals of INDRA 



The total light output friendly formula (9) was used to perform a compar- 
ative study of the CsI(Tl) crystals of INDRA. At forward angles (rings 2 - 
9), the experimental values were obtained from elastically scattered light ions 
on various targets (C, Al, Au). The light ions were produced by secondary 
beams with atomic number Z = 1 - 6 and mass number A precisely identi- 
fied. The threshold energy per nucleon required to generate 5 - rays which 
contribute to the scintillation is an intrinsic characteristic of the Csl crystals. 
Consequently, after a grid search over all the modules of the 8 forward rings, 
the parameters 03, and hence too, were fixed. The crystals belonging to one 
ring have the same size, the same shape and similar Tl (as well as eventually 
defect) concentrations. The gains of the associated PMT are the same. For 
this reason, one would expect that the parameters ai and 02: 

. . Aac-^V 1 
«1 OC /geom X /pMT X — 02 CX 



AAeA^A + AoA^o' AAeA^A + AoATo' 



take nearly the same values for the modules of one ring; the quantities AacA^Aj 
AdA'd were defined in the previous paper [3] and /geom, /pmt are factors 
connected to the light collection (geometry of the crystal) and to the PMT 
gain, respectively. Actually, this is the situation, as shown in Fig. 7a), b) for 
the modules of ring 2. 

If AdA^d/AacAa 1, the parameter ai//geom is niainly related to the associ- 
ated PMT gain. The geometrical light collection factor /geom is proportional 
to the response of the crystal at ^^^Cs source 7 - ray irradiation measured 
with the same PMT for all the crystals of INDRA. Averaged over the mod- 
ules of the same ring, the parameter < ai > / < /geom > plotted versus 
the ring number, follows the approximately known values of the PMT gain 
(provided by the manufacturer) as shown in Fig. 7c). The correlation of the 
two parameters < Oi > and < 02 > may be followed in Fig. 7d), if < Oi > 
is previously corrected for the geometric and gain factors. Obviously, there 
is no mathematical correlation: two groups of detectors are put in evidence. 
They may correspond to the different concentrations of the Tl activator - from 
200 ppm to 2000 ppni - (and eventually of the other crystal imperfection) 
which could appear during the CsI(Tl) crystal growth through the Bridgman 
method [20]. At backward angles of INDRA (rings 10 - 17), only light charged 
particles, available from secondary beams at GANIL, could be used in the fit 
procedure to get the results presented in Fig. 7c), d). For unitarity, in these 
last two plots, we have also restricted the forward angle employed data to light 
charged particles only. 
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Fig. 7. a), b) The fit parameters ai, 02 used in the simple formula (9) against the 
module number for ring 2 of INDRA. c) The parameter < ai >, averaged over 
all the modules of each ring and corrected for the light collection, versus the ring 
number; it follows the PMT gains - hatched area, d) The quenching parameter 
< 02 >, averaged over the modules of each ring, versus the parameter < ai > 
corrected for light collection and PMT gain, averaged over the modules of the same 
ring. The error bars arc smaller than the symbols because only statistical errors on 
the parameters have been considered. 



The above statistics performed over more than 300 CsI(Tl) crystals belonging 
to INDRA array and the related considerations enforce the consistancy of the 
RNQM and of the exact or approximate light output expressions obtained in 
its framework. The parameters 03, 04 and, to a certain extent, 02 may be con- 
sidered as intrinsic characteristics of the Csl(Tl) scintillators used in nuclear 
physics applications. Insofar, the rather simple formula (9) of the total light 
output, together with the recommended values in Table lb) constitute a pow- 
erful tool to be used in heavy ion experiments. Starting from easily to obtain 
light charged particle experimental data (excluding protons if their total light 
signal is not directly measured), they allow to accomplish useful applications 
such as heavy fragment identification and their energy determination, in case 
of lack of directly related information. An example is given in the following 
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subsection. 



4-4 Reaction product identification and energy calibration of the backward 
angle modules of INDRA 

There are two major difficulties when the CsI(Tl) scintillators of INDRA 
placed at polar angle above 45° are exploited. Firstly, in a ^Eic — Qo map the 
identification is perfect for low atomic number Z of the fragments {Z < 15). 
For higher Z values the identification of the atomic number becomes uncer- 
tain, because of the poor statistics and the lower IC energy resolution, as 
compared to the Si detectors. Secondly, at large angles it is quite difficult to 
obtain data of elastic scattering for fragments. Insofar, the calibration of the 
module belonging to the same ring is performed in two stages. 

In the first stage, we consider the CsI(Tl) partially obturated by the cali- 
bration telescope CT (Si(80 fim) - Si(Li) 2 mm), one per ring. Most of the 
fragments {Z >4) are stopped in the Si(Li) and therefore identified in the Si - 
Si(Li) telescope. The corresponding energy spectrum is built for each atomic 
number Z. There are energetic hght charged particles (> 20 AMeV) passing 
through the silicon telescope and stopped in the CsI(Tl) crystal coupled be- 
hind. The calculated [9,10] residual energies are the "true" energies deposited 
in the scintillator. Relied on the analytical expression (9) of the light output, 
and using the recommended values from Table lb) for parameters 02,03 and 
04, the parameter oi is determined by a fit procedure. By means of these 
parameters, the energy spectrum in the scintillator is built for each of the 
light charged particles punching through the silicon telescope. Put together 
the silicon telescope and CsI(Tl) scintillator spectra provide the whole energy 
spectrum for a given particle. These spectra and those of fragments up to 
Z—15 stopped in the Si - Si(Li) telescope are the reference spectra for the re- 
spective ring, concerning an energy range: £"111111; -E-max for each type of reaction 
product. At the same time, the parameters allow to perform Z identification 
in a AEic — Qo map for the whole range of atomic number of light charged 
particles and fragments passing besides the silicon telescope and entering the 
CsI(Tl) coupled behind. The energy spectra of these reaction products were 
compared with the reference spectra; the superposition is quite good. 

In the second stage all the other CsI(Tl) scintillators of the ring are calibrated. 
The energy spectra of a given Z have to be identical to the corresponding refer- 
ence one (the trigger used in the experiments did not break the axial symmetry 
of INDRA). An energy spectrum, starting from the light response spectrum, 
is obtained by stretching it between -Emim -^max in order to reproduce the ref- 
erence spectrum. A minimizing procedure based on the MINUIT package 
from CERN library is used. It provides the parameters Oi, for the considered 
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Fig. 8. Energy spectra of light charged particles and fragments in one of the CsI(Tl) 
crystal of ring 10 of INDRA (symbols). They are compared to the same charged 
reaction product energy spectra of the reference module (lines) of this ring. 

module. The parameters Oi (i=l,4) allow afterwards a complete Z identifica- 
tion in the AEic — Qo map. The energy spectra of the fragments stopped in 
the respective scintillator arc determined too [23]. The good quality of the 
procedure is illustrated in Fig. 8, where energy spectra of different particles 
stopped in one of the scintillators of ring 10 of INDRA are compared to the 
reference spectra of the same ring. 



4-5 And if the 5 - rays would be neglected? 



We have to stress once again the importance of taking into account the 5 - 
ray effect in the light output, especially for reaction products heavier than 
intermediate mass fragments {Z > 15). By neglecting it, only the first term 
in the right side of Eq. (9) would appear. The results remain reasonable for 
Z < 15 but with values of the fit parameters which change significantly as 
compared to the case where the knock-on electrons were taken into account. 
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The quality of the fit is drastically degraded for fragments heavier than Z—15. 
So are particle identification (Fig. 4) and energy calibration. 



5 Conclusions 

The data obtained with the INDRA array (large range in Z and £"0) provided 
a good opportunity to proceed to a more basic study of the light output 
of CsI(Tl) scintillators and to derive suitable calibration and identification 
procedures. Starting from the fast and the slow parts of the hght output 
(fractions of the total light response integrated in the corresponding time 
gates) , it was possible to rebuild the integral of the signal. 

Under suitable approximations, the expression of the total light output de- 
rived in the preceding paper [3] may be analytically integrated. Even if up to 
3% of the accuracy may be lost, the fact presents the huge advantage of ex- 
tremely short computing time. The derived expressions, easily to handle, were 
successfully applied for fragment identification in Ai?si,ic ~ Qo maps and for 
the energy calibration of the CsI(Tl) scintillators. At forward angles, where a 
Si detection layer exists, these applications lead to an important reduction of 
the computing time. At backward angles, where two problems exist: fragment 
identification and energy calibration of the CsI(Tl) crystals, the above proce- 
dure plays an even more important role and corresponds to the optimum way 
we have found to solve these two tasks. 

A comparative study of the CsI(Tl) scintillators of INDRA has shown that 
the model parameters are meaningful quantities, related to the light collection 
and the PMT gain, to the activator and eventually crystal imperfection con- 
centrations and to the 6 - ray production energy threshold. Except the gain 
parameter, all the others are characteristics of the usual CsI(Tl) scintillators. 
Their averages, performed over the 324 CsI(Tl) crystals of INDRA, allowed to 
find reliable, recommended parameter values. Together with the related total 
light expression, they constitute good implements for energy calibration and 
heavy ion identification applications in heavy ion physics experiments. 
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